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We study the superconducting instabilities ol a single species ol two-dimensional Rashba-Dirac 
fermions, as it pertains to the surface of a three-dimensional time-reversal symmetric topological 
band insulator. We also discuss the similarities as well as the differences between this problem and 
that of superconductivity in two-dimensional time-reversal symmetric noncentrosymmetric materials 
with spin-orbit interactions. The superconducting order parameter has both s-wave and p-wave 
components, even when the superconducting pair potential only transfers either pure singlet or 
pure triplet pairs of electrons in and out of the condensate, a corollary to the nonconservation of 
spin due to the spin-orbit coupling. We identify one single superconducting regime in the case of 
superconductivity in the topological surface states (Rashba-Dirac limit), irrespective of the relative 
strength between singlet and triplet pair potentials. In contrast, in the Fermi limit relevant to 
the noncentrosymmetric materials we find two regimes depending on the value of the chemical 
potential and the relative strength between singlet and triplet potentials. We construct explicitly 
the Majorana bound states in these regimes. In the single regime for the case of the Rashba-Dirac 
limit, there exists one and only one Majorana fermion bound to the core of an isolated vortex. In the 
Fermi limit, there are always an even number (0 or 2 depending on the regime) of Majorana fermions 
bound to the core of an isolated vortex. In all cases, the vorticity required to bind Majorana fermions 
is quantized in units of the flux quantum, in contrast to the half flux in the case of two-dimensional 
Pj. ± ipy superconductors that break time-reversal symmetry. 
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I. INTRODUCTION 

Bi2Sc3 is an inversion-symmetric layered Ipaijd insula- 
tor with a bulk gap estimated to be 0.35 eV.Eru Density- 
functional theory predicts that Bi2Se3 supports a single 
Rashba-Dirac cone of gapless surface states, a prediction 
that has been veiiBed using angle-resolved photoemis- 
sion spectroscopy.EfQ This remarkable attribute of Bi2Sc3, 
which has otherwise onlyJpeen observed in the insulating 
alloys Bii_2;Sb^ so farJj'El is the defining property of a 
three-dimensional (3D) tinjMvrevcrsal symmetric (TRS) 
topDlogical band insulator. lHj In a recent work, Hor et 
aZ.E3have reported the observation of strongly type II su- 
perconductivity in Cu2:Bi2Se3 below 3.8 K when Cu is 
intercalated between the Bi2Se3 layers. They have also 
proposed to use Cua;Bi2Se3 as a mean to induce super- 
conducting correlations for the TRS topological surface 
states by the proximity effect. 

The surface states in a 31? TRS topological band in- 
sulator are reminiscent of the Bloch states of graphene 
in that, in both cases, their density of states vanishes 
linearly at the so-called Rashba-Dirac point.EJ However, 
they differ in a fundamental way from those in graphene. 
For example, the surface of Bi2Se3 supports one Rashba- 
Dirac cone as opposed to two in graphene. This difference 
is a manifestation of the fact that inversion symmetry is 
maximally broken on the surface of Bi2Se3 in that the ki- 
netic energy is dominantly of the Rashba type, whereas 
the spin-orbit coupling is for all intent and purposes neg- 
ligible for graphene. Consequently, the surface states of 



a 3D TRS topological, hand insulator are not localized by 
weak TRS disorder ii^llS whereas Anderson localization 
rules in graphene. 

Another difference with graphene, as we shall show in 
this paper as a warm up, is that all states in the Rashba- 
Dirac sea contribute to the Pauli magnetic susceptibil- 
ity, which is anisotropic in that the in-plane and out-of- 
plane components differ by a factor of 2. For compari- 
son, the Pauli magnetic susceptibility is isotropic in spin 
space and proportional to the density of states at the 
Fermi surface in any electron gas (including graphene) 
with small breaking of the spin-rotation symmetry (SRS). 
This anisotropy and the fact that the Pauli susceptibil- 
ity does not only depend on the density of states at the 
Fermi level could potentially be used as a simple diagnos- 
tic of a limit in which the Rashba coupling is the largest 
energy scale. 

The main emphasis of this paper will be on the super- 
conducting instabilities of the surface states in a 3D TRS 
topological band insulator and on those in close relatives, 
i.e., two-dimensional (2D) TRS noncentrosymmetric su- 
perconductors in a regime that has been little studied so 
far. The theoretical studies of noncentrosymmetric su- 
perconductors with TRS usually assume the hierarchy of 
energy scales 



t > a > A 



(1.1) 



where t is the inversion-symmetric band width, a > 
is the spin-orbit coupling that preserves TRS but breaks 
SRS, and A > is the single-particle superconducting 
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The regime 

a>A>t (1.2) 

is the one that apphes to intrinsic superconducting insta- 
bihties of the surface states in a 3D TRS topological band 
insulator. We will address the question of whether inter- 
esting phenomena associated to s upe rcond uctiv ity occur 
upon exchanging the hierarchies ( |1 . l| ) and ( |l.2| ). 

In the same way that (2D) TRS band insula- 
tors have -been, classified according to their topological 
character Bogoliubov-de-Gennes (BdG) supercon- 
ductors have also been given topological attributes when- 
ever they support gapless boundary states in confined 
geometries.Eil^Ej A necessary (but not sufficient) condi- 
tion for a 2D TRS superconductor to be topologically 
non-trivial is that it is noncentrosymmetric, or, equiva- 
lently, that it breaks SRS. According to Refs. ^ and |3^, 
a sufficient condition is that, for any weak and local TRS 
static disorder, a 2D TRS superconductor in an infinitely 
long strip geometry supports an odd number of Kramers' 
doublets of gapless edge states of which at least one wave 
function is extended along the edge (see also Refs. p^]3^ 
for varying alternative criteria). 

Applying this definition of a topological superconduct- 
ing phase to the superconducting instabilities of surface 
states in a 3D TRS topological band insulator immedi- 
ately leads to a paradox: What is the meaning of the 
boundary of a boundary? A more meaningful question 
to ask might be: What are the spectral properties of 
TRS-breaking vortices if the surface states in a 3D TRS 
topological band insulator support a type II supercon- 
ducting order? Do they bind mid gap states generically, 
zero modes in particular, or not? These are questions 
that we address in this paper. 

Defects in a type II supercanductor are vortices. 
On the one hand, Caroli et a/£3have shown that vor- 
tices in an s-wave TRS and SRS superconductor bind 
nonvanishing-energy bound states with a mean level 
spacing of the order of the superconducting gap squared 
divided by the Fermi energy. On the other hand, Jackiw 
and Rossi in Ref. ^ found a single bound state that is ex- 
ponentially localized around the core of a unit-flux vortex 
in a 2D s-wave relativistic superconductor with a vanish- 
ing density of states (Rashba-Dirac point). The energy of 
this bound state is precisely pinned to the Fermi energy 
(see also Re f. |39| for the corresponding index theorem and 
Refs. ^ and |4l|for examples of nonrelativistic zero modes 
bound to vortices) . A midgap state bound to the core of 
a vortex docs not carry an electric charge, for it is an 
eigenstate of the generator of the particle-hole symmetry 
(PHS) obeyed by any BdG Hamiltonian. It is thus charge 
neutral and as such is a physical realization of a Majo- 
rana fermion. Majorana fermions were also found to be 
exponentially localized to the core of a vortex in ap^±ipy 
type II superconductor by Read and Green and by Ivanov 
in Refs. ^ and respectively. More importantly, they 
showed that these Majorana fermions obey non-Abclian 
braiding statistics. Theoretical proposals to nucleate Ma- 



jorana fermions have been made relKipg on 2D TRS non- 
centrosymmetric superconductors,CjE3 or on proximity 
effects at the 2D interface betweea. hand insulators, su- 
perconductors, and ferromagnets.cZrE3 

We will show in this paper that, when the dispersion is 
Rashba-Dirac like, there is a single zero mode bound to 
the core of an isolated vortex with unit circulation, and 
thus a single Majorana bound state. The mechanism, 
in the case-of singlet pairing, is precisely that of Jackiw 
and Rossi.Ea This zero mode remains for arbitrary ratios 
of triplet and singlet pairing, with the pairing potentials 
Aj and Ag , respectively, and also as the chemical poten- 
tial ^ is varied. The stability of a singly degenerate zero 
mode is guaranteed in a system with particle-hole sym- 
metry in which the zero mode is isolated from the contin- 
uous spectrum by a finite energy gap. Therefore, study- 
ing gap-closing surfaces in the parameter space of the 
coupling constants characterizing the theory is of crucial 
importance in identifying the stability of the Majorana 
modes as well as the phase boundaries between different 
topological phases. 

In this paper, we compute the conditions for the closing 
of the gap in A^/A^-fj, space by exploring a one-to-one 
mapping to the normal-state dispersion relation, in which 
a function of the ratio A^. / A^ serves as a reparameteriza- 
tion of the magnitude of the momenta in the dispersion 
relation. Thereto, we show that there are as many lines 
in At/Ag-/i space at which the gap closes as there are 
branches in the dispersion relation. But in the case of the 
Rashba-Dirac dispersion, it is possible to go from one side 
of a gap vanishing line to another without crossing it by 
going through the point at infinity (Ag — 0). Thus, there 
are not two distinct phases separated by a transition in 
this case, but there is a single phase instead. 

In the Fermi limit relevant to 2D TRS noncentrosym- 
metric superconductors, we find that the conditions for 
the closing of the spectral gap do separate two gapful 
phases. These two regimes are those in which either the 
singlet or the triplet pairing controls the physics. The 
detailed shape of the phase boundaries is dictated by 
the normal-state dispersion relation. The presence of 
the TRS spin-orbit coupling leads to interesting effects 
at certain values of the chemical potential, for example, 
re-entrance to the phase dominated by singlet-pairing 
physics even when A^/Ag is large. We find two Ma- 
jorana zero modes bound to an isolated vortex in the 
triplet controlled phase but they disappear in the singlet 
controlled phase. We find that the vortices that bind this 
pair of Majorana zero modes have unit flux, as opposed 
to the half- vortices needed in the case of two-dimensional 

± ipy superconductors that break time-reversal sym- 
metry. The physical reason for this difference is that, 
when TRS holds, the spin-resolved pairing amplitudes 
A^l and A^^ are not independent, and thus one cannot 
introduce vorticity in one but not the other, which is 
the case for the half vortices in the p^ ± ipy supercon- 
ductors. Therefore the pair of Majorana fermions that 
we find for full vortices in the triplet case is distinct from 



those found by Read and GreenB and Ivanov.@ The pair 
of Majorana fermions that we find is not robust to a 
generic weak perturbation that breaks translation invari- 
ance, for these Majorana fermions are not related to each 
other by the operation of time reversal. This pair of Ma- 
jorana fermions is thus unrelated to the one introduced 
by Qi et al. in Rcf. ^ as a mean to identify the triplet 
dominated TRS phase as a nontrivial 2D topologi- 
cal superconducting phase. We conclude that, although 
both superconducting regimes can be distinguished by 
the even number of Majorana fermions that an isolated 
TRS-breaking vortex binds, this distinction is not topo- 
logical, for it is not robust to static disorder, for example. 

In addition to this interesting interplay between the 
singlet and triplet pair potentials for the existence of Ma- 
jorana fermions, superconductivity for surface states in a 
3D TRS topological band insulator and noncentrosym- 
mctric materials has other curious properties, one of 
which is the following. Because spin is not a good quan- 
tum number, even if the pair potential contains, say, only 
the singlet component, the condensate will nevertheless 
have triplet correlations. For example, the pairing corre- 
lation for electrons (c^-|.c_j.-f-) 7^ (and (c;^.^c_^^) 7^ as 
well) for /i 7^ 0, even if = 0. One measurable conse- 
quence is that, generically, these superconducting states 
would lead to detectable Josephson currents when con- 
nected to either conventional s-wave or p-wave supercon- 
ductors, manifesting the fact that they have both types 
of correlations (even if only the singlet pair potential 
is nonzero). 

The paper is organized as follows. We define the model 
in Sec. We show in Sec. ^ that the normal-state Pauli 
magnetic susceptibility has the remarkable property that 
it depends on all states in the Fermi sea, not only on 
the density of states at the Fermi level, and that it is 
anisotropic for the surface states of a 3D TRS topolog- 
ical band insulator. The dynamical Pauli susceptibility 
tensor in both the normal and superconducting states 
encodes rich magnetoelectric effects that are responsi- 
ble for the spin-Hall effect among others. We study in 
Sec. 1^ sclf-consistently the interplay between the sin- 
glet and triplet components to the superconducting in- 
stabilities of the surface states in a 3D TRS topological 
band insulator or in 2D noncentrosymmetric materials. 
The generic mean-field phase diagram for a TRS two- 
band BdG Hamiltonian in the isotropic continuum limit 
that interpolates between the regimes (1.1) and ( |1.2D is 
constructed in Sec. 0. We find that the Rashba-Dirac 
limit t/a = 0, that pertains to the surface states of a 
3D TRS topological band insulator, is singular in that 
there ex ists o nly one single phase in the phase diagram. 
In Sec. VIA, we construct explicitly the single Majo- 



rana fermion bound to the core of an isolated vortex in 
the superconducting phase of the surface states of a 3D 
TRS topological band insulator. In Sec. VI B , we con- 
struct explicitly the pair of Majorana fermions bound to 
the core of an isolated vortex in the triplet dominated 
superconducting phase of a 2D noncentrosymmetric su- 




FIG. 1: (Color online) (a) Schematic picture of the surface 
states of the topological insulator Bi2Se3 (black lines). The 
chemical potential is far from the Rashba-Dirac (nodal) 
point and close to the conduction-band continuum (upper 
gray region), while the Rashba-Dirac (nodal) point is close 
to the valence-band continuum (lower gray region), (b) Left: 
one-dimensional cut o f th e dispersion of the Rashba-Dirac 
model defined by Eq. (2.1). In particular, we study the case 
where ^ is close to the Rashba-Dirac (nodal) point rather 
than close to the energy cutofi' ±A that defines the onset of 
the conduction band and the valence band. Right: the expec- 
tation values of the electron spins are perpendicular to their 
momenta and oriented i n opp osite direc tions for the upper 
and lower cone [see Eqs. (3.3b) and (3.3c) at 5 = 0]. 



perconductor. We conclude with Sec. VII. 



II. 



DEFINITION 



Normal state 



In the continuum limit and in the single-particle ap- 
proximation, we define the single-node Rashba-Dirac 
Hamiltonian 



nisur 



hv 



RD 



|fc|<A. (2.1) 



The Rashba-Dirac velocity is Wj^q and we restrict the 
momentum hk by the cutoff A/wj^j-, beyond which the 
surface states of a TRS topological band insulator merge 
into the bulk states. The two-dimensional momentum 
hk = h{k-^^^ couples to the Pauli matrices cr = (ctj^, (T2). 
These Pauli matrices act on the internal space of the spin- 
1/2 degrees of freedom carried by the surface electron 
(hole) in the laboratory frame of reference. This coupling 
between the electron (hole) crystal wave vector and the 
spin of the electron (hole) prevents conservation of the 
electron (hole) spin. However, TRS is con serv ed so that 
the linear dispersion that follows from Eq. (2.1) is twofold 
Kramers degenerate. 

For the surface states of Bi2Se3, the RashbaJDirac ve- 
locity is measured to be « 5.1 x 10^ m s~^.cl Further- 
more, the Rashba-Dirac energy (the energy measured at 
the Rashba-Dirac point) e^j-, « 0.3 eV is close to the in- 
sulating band gap of 0.35 eV for the bulk states in which 
the surface states mergcQ Hence, the superconducting 
gap A w 3 X 10~^ eV in intercalated Cua:Bi2Se3 is minute 
compared to in Bi2Sc3. If there are Rashba-Dirac 
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surface states in Cu2;Bi2Sc3 involved in pairing correla- 
tions or if there are Rashba-Dirac surface states in Bi2Se3 
involved in pairing correlations induced by the proximity 
to superconducting Cua;Bi2Se3, they are likely to be far 
away from the Rashba-Dirac point. On the other hand, 
we take the point of view that it is only a matter of time 
before a way is found to tune the chemical potential of the 
TRS topological surface states through the Rashba-Dirac 
point (substituted CuyBi2_ySe3 might be a candidate). 
Hence, one goal of this paper is to characterize pairing 
correlations among the surface states of a TRS topologi- 
cal band insulator upon tuning the Fermi energy through 
the Rashba-Dirac point (see Fig. |l|). 



B. Mean- field superconducting state 



We shall compare our study of Eq. (2.1) with that of 
the two-dimensional Rashba tight-binding model 



-1/2D 



9k - 



k e BZ. 



(2.2a) 



Here, ctq is the unit 2x2 matrix in spin space and the 
wave vector k is restricted to the first Brillouin zone (BZ). 
It describes the hopping on a square lattice with the SRS 
dispersion 



~ —2t (cos ki + cos fca) , t > 0, 
and with the Rashba spin-orbit coupling 



9k 



a 



— sin h 
sin ki 



(2.2b) 



(2.2c) 



say. Our convention throughout this paper will be that 
a > 0. 

An important difference between the surface states of 
a 3D TRS topological band insulator and the Rashba 
tight-binding states is that the surface states span an 
odd number less of Fermi surfaces. This is a manifesta- 
tion of the fermion doubling that occurs when attempting 
to regularize a ZJ-dimensional single Rashba-Dirac cone 
by a D-dimensional tight-binding model. The fermion 
doubling can only be eliminated by the addition of the 
Wilson term 

■Hjf^ := % (2 - cosfc^ - cosfcj) fs, tyj :$> t + a, (2.3) 



to the tight-binding Hamiltonian (2.2a) at the cost of 
breaking TRS. 

Hamiltonian (2.1) is scale invariant. It then follows 
that the density of states per unit area is proportional to 
the absolute value of the chemical potential ^ and van- 
ishes at the Rashba-Dirac point = 0. The effects of 
this scale invariance on charge transport, including the 
orbital effects of a magnetic field, are identical to those in 
graphcnc in the single Rashba-Dirac cone approximation, 
if the Zceman coupling to a magnetic field is ignored. In 
Sec. [II, we are going to study the inherently strong ef- 
fects of the spin-orbit coupling on the Pauli magnetic 
susceptibility. However, before doing so, we want to in- 
clude the possibility of a superconducting instability that 
we first treat at the mean-field level without imposing the 
condition of self-consistency. 



We rewrite the c onti nuum Hamiltonian ( |2.l[ ) or the 
lattice Hamiltonian ( |2.2| ) in the language of second quan- 
tization. For simplicity, we choose a tight-binding nota- 
tion. Reverting notation to the continuum is straight- 
forward. We thus introduce the spinor ■0^ — ("-It' '^kij 
for electrons in the spin basis of the laboratory frame of 
reference and the spinor ^/i^ = ^a^^, '^I-) helicity 
basis defined below. This gives 



feeBZ fceBZ 
'^a-.k = - m) Ctq + Sfc • cr. 



(2.4a) 



Hok 



k+ 







Cfc. 

The single-particle dispersion is here given by0 

4± -£fc-M±l9fcl, (2.4b) 

while the transformation between the laboratory basis 
and the helicity basis is given by the unitary 2x2 matrix 



1 / 1 



1 



1/1 1 

—/= 3fci+'gfe2 gfei+'gfc2 



(2.4c) 



whereby 

- ^i^l V-fc - HfeV^fc, Ho;k = ^ino;k^k- (2.4d) 
The (mean-field) BdG Hamiltonian is defined by 



keBZ V ^k "O;-* 



f '^o-.k ^k , ^ ^ ^ 



-^o.k ^k 



keBZ 



where the bispinors 4*^, and arc given by 



respectively. We have chosen to construct the bispinors 
$^ and $^ from the spinors V'fe and ^/j^, respectively, and 
their time-reversed partners. Thereby, we have to take 
care of the action of the time-reversal operation T on the 
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laboratory and the helicity single-particle states labeled 
by the wave vector k and the indices s =t, i and A = ±1, 
respectively. For the laboratory basis, it is 

r|fct> = + |-fc;), r|fe;> = -| -fet), (2.6) 

i.e., it is off-diagonal in the laboratory spin basis. For 
the helicity basis. 



|fc+> = 
|fc-) = 
together with 



1 

1 



(lfct> 
(lfct> 



fc; )) 
fc; )) 



imply that it is 

TlkX 



Ae-''^-'= I - fcA ; 



(2.7a) 



(2.7b) 



(2.7c) 



i.e., it is diagonal in the helicity internal space but with 
the wave vector and helicity-dependent eigenvalue Ae"^*" 
that is odd under the inversion k —k. Hence, the 
bispinors $^ and follow. 

We parameterize the 2x2 pair-potential matrix by 



A, 



(2.8a) 



in the laboratory frame for the spin degrees of freedom. 
PHS, which embodies Fermi statistics within the BdG 
formulation, demands that it is an antisymmetric opera- 
tor, i.e.. 



TRS imposes the conditions 



-fc: 



d-k = -dlk- 



(2.8b) 



(2.8c) 



Throughout this paper, we consider Cooper pairs 
made of time-reversed helicity single-particle states from 
Eq. (2.7c). Hence, we take the 2x2 pair-potential matrix 



A, 



A 



A, 



(2.9a) 



to be diagonal in the helicity basis, and it then follows 
that 



A — A* 



A — A* 



(2.9b) 



as a consequ ence of TRS. Furthermore, we find with the 
help of Eq. (2.5b) the 4x4 Hermitian matrix (the com- 
plex notation z = x + iy and z = a; — iy is occasionally 
used) 



/^^fc - 9k 

9k £fc-M ^ 

c.c. c.c. 

c.c. c.c. 



t,k 



s,fc 



A 



S,fc 

9k 



-9k 



(2.10a) 



where we recall that ip/^ := arg gj^ and 



A. , = i ( A 



^k+ 



A 



= A* 



At_fc = -(A^^-A 



dk= Afc+ ^ Afc- 



= A, 



9k 

\9k\ 



t.-fcj 



= -d* 



(2.10b) 



in the representation of Eq. ( 2.5a ). 

The fact that the vector is parallel to is a 
consequence of our assumption that Cooper pairs are 
mad e of t ime-reversed helicity single-particle states, i.e., 
Eq. ( ^.9a ) . This assumption is justified if the pairing in- 
teraction preserves the symmetry of the noninteracting 
Hamiltonian. Following|-tlic literature on noncentrosym- 
metric superconductorsJH3 wc arc thus assuming that the 
symmetry of the noninteracting Hamiltonian is preserved 
by the self-consistent inclusion of the pairing interaction. 

We also demand that Hamiltonian ( ^.lOa ) is single val- 
ued in the BZ. This restricts the triplet pairing A^ ^ to 
vanish at least as fast as 



Wk\ 

hm - 

l9j->0 \9k\ 



< c 



(2.10c) 



for some number c larger than or equal to 0. With our 
choice of gauge, A^._|_ and A^_ or. cquivalently, A^ ^ and 
Aj. J. arc real valued. In Sec. where we study TRS- 
breaking vortices, we revert instead to complex order pa- 
rameters to accommodate twists in the phases of the sin- 
glet and triplet pair potentials. Finally, we observe that 
the pair-potential eigenvalues 



AfcA = \k 



9k 

\9k\ 



A = ±, 



(2.11) 



transform according to the same irreducible representa- 
tion of the space group. For example, in the isotropic 
continuum limit with s-wave pairing they are functions 

of |fe| only. 

The BdG Hamiltonian ( |2.5a| ) is of the form 

A=± A=±fceBZ 

Hk\ = ^kx^lxa^x + AAfc;, (e"^-'=a_fcAafcA + H.c.) . 

(2.12) 

The mean- field ground state is the state |Tj^f) anni- 
hilated by H. It is obtained as the direct product 



^mf)' where |T^f) is annihilated by 



for each of the helicities A 



Each helicity supports quasiparticles obeying the PHS 
(relative to the chemical potential) dispersion iE^,^ with 



^kx - 



^kX 



A^ 

^fcA' 



A = ±. 



(2.13) 



However, the ground states |T^f) and are not inde- 

pendent since they arc tied to each other by TRS. Indeed, 



6 



TRS implies that the relative phase of the pairing poten- 
tials Aj^.^ with the helicities A = ± is locked to be or 



TT, as follows from the transformation law (2.7c), i.e., 

(2.14) 

To construct ITj^^f), we perform a Bogoliubov trans- 
formation for each helicity index A = ± independently. 
Thus, for each helicity A = ±, we define 

lk\ ■= Ukx flfcA - Vkx alfcA (2.15a) 
with the complex- valued coefficients ?7^_|^ and V/,^, 



fcAl 



-2" 



fcA 



fcAl 



i 

= o ( 1 - ) ' (2-15b) 



fcA 



El 



fcA 



kX 



'kX ^kX 

Under this transformation 



fcA 



H = J2 ^kX llxlkX 



(2.16) 



fceBZ A=± 

The mean-field ground state is then 

-.t „t 

\ykX ''fcA' 

A=± k 

provided 

7fcA|T,„t) = 



T„,f> = n n {^kx + V^A^A^lfcA) |0) (2-17a) 



(2.17b) 



holds for all k and all A = ±. 

By construction, the mean-field ground state ( 2.17 ) is 
TRS. SRS is, however, broken. Consequently, 



C^mf |C-fc-|-Cfct I "^mf ) 



C^mf |C-fc4.Cfc4_|Tj,jf) 
C^mf I C-fc-|-Cj.j_ I T„jf ) 




(2.18) 



are generically nonvanishing (one exception is the 
Rashba-Dirac limit e^^ = at the Rashba-Dirac point 
11 = 0) even though the pair potential may be purely 
singlet when 



A, 



A 



or purely triplet when 
A, 



-A 



fc-- 



In a superconducting state that preserves SRS, the 
ground state has no spin correlations other than that 
of the pair condensate. 



III. 



SUSCEPTIBILITY 



A. Static and uniform Pauli magnetic 
susceptibility at T = in the normal state 

We are after the Pauli magnetization per electron in- 
duced by the Zeeman coupling oc —BiOi — i?2'''2 ^ B^a^, 
where it is understood that Cg is the third Pauli matrix 
and the in-plane components of the magnetic field are 
Bi and B2 while the out-of-plane component B^ is taken 
along the spin quantization axis in the laboratory frame 
of reference. 

To obtain the Pauli ma gnet ization per electron at 
T = 0, we start from Eq. (2J) with the Zeeman cou- 
pling added 



o/sur 
i^Bk 



(^%D^1 - ^1 " {^^WK + -Si) - BgCTg, 



(3.1) 

compute the expectation value of the spin operator 
ha-^ 2 3/2 for all the Bloch states, and sum these expec- 
tation values up to the chemical potential [i. The Pauli 
susceptibility per electron then follows by differentiation 
with respect to B^ ^ 3 followed by setting i?^ 2 3 = 0. We 
set h = Upjj-, = 1 to simplify notation. 
As long a& B\ + B\ > 0, the eigenvalue 



e±(fc) = -M ± V (^-1 - ^2) + (^2 + Si) + Bl (3.2a) 
has the eigenstate 



1 



AA±(fc) 
with the normalization 



— ifc^ — + iSj 

e±(fe)+M + S3 



(3.2b) 



AA±(fc) := ^2[e±(fc) + /x] [^±(fe)+;, + B3]. (3.2c) 

We observe that the effect of in-plane magnetic fields is 
to translate the Fermi sea. The spin expectation values 
in the Bloch states are 



B, 



Bi 



B 



T^AF{B/-k,B/k) 



(3.3a) 



(2.19) ^^>l{k)a^^{k) 



(2.20) 



2{k + B) (e±(fc)-t-/.) 



S|2 

B 



k + B\ 



B 



[F {B/k,B/k) 



(3.3b) 
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and 
^']t(A;)a-«'±(fc) 



2{k + B) (C±(fc)+/^) 



B,=0 



= T 



|fc + B|2 
k + B 



(e±(fe) 



|fc + B 



Ifcl 



F{B/k,B/k) 



Here, we have introduced the complex notations 



'1 T 10'2' 

for the Pauli matrices, 



i^i, 



for the momenta, and 



B^B, 



ifi 



2i 



(3.3c) 



(3.4a) 



(3.4b) 



(3.4c) 



for the in-plane components of the magnetic field. We 
have also introduced the real-valued function 



F{z,z) 



1 



1 



{z + z) 



(3.4d) 



that comes about to first order in an expansion in pow- 
ers of the components of the magnetic field. The mag- 
netization per electron is obtained by integrating over 
all single-particle energies up to the chemical potential 
jjL. We conclude that the Pauli magnetic susceptibility 
tensor per electron is 



Xab « X 



7r(A-|/^|), if a =1,2, 



27r(A-|//|), if a = 3, 



(3.5) 



in the noninteracting approximation and at T = 0. 

The Pauli magnetic susceptibility per electron 
also holds for the Rashba tight-binding Hamiltonian 
with minor modifications provided the limit </q — > 
taken, {N is the number of lattice sites) 



3^) 
2.2) 



is 



f 1/2 ^ 

N 

M<\9k\ 



1 

Is J' 



Xab = (^ab X S 



lMl<lafcl 



if a = 1,2, 



if a = 3. 



(3.6) 



In the opposite limit a/t — > 0, we cannot use the lat- 



tice counterpart to Eq. (3.3a) to compute X33, since our 
choice for the spinor representation is singular in this 
limit. We can however use the lattice counterparts to 
Eqs. ( 3.3b ) and ( [3.3cD to compute Xii and X22- By 
isotropy, we then recover the conventional Pauli magnetic 
susceptibility 



Xab « 5^bVp{^x) 



(3.7) 



where v-p^fi) is the density of states per electron and per 
spin of the dispersion e^. This result remains true to first 
order in a/t. 



B. Dynamical Pauli susceptibility in the 
superconducting state 

Another remarkable consequence of the spin-orbit cou- 
pling is that charge-density and spin-density fiuctuations 
are coupled., bot^i in the normal and in the superconduct- 
ing state.lnl£irE2U|Clie spin-Hall effect is a consequence 
of this coupling .OH To quantify this statement, we in- 
troduce the susceptibility tensor in the superconducting 
state 



(Xoo), 

{X0d)q 
{Xbo)q 
{Xbd)q 



— T 



tr [G, 



0;fc^03^0;fc+g^03] ' 



' [^0;k^boGo-k+q^do] > 



(3.8a) 



where the indices b and d run over the values 1,2, and 3, 
and 



^,1^ = 0, 1,2,3, 



(3.8b) 



with the unit 2x2 matrix Tq and the Pauli matrices 
T acting on the particle-hole two-dimensional subspace. 
According to t he me an-field Hamiltonian in the supercon- 
ducting state ( 2.5a ), the single-particle Green's function 
is 



C!o;k [^i^n-'^OO + i^k ^ M) ^03 + 9kl^l3 

+ ^s,fe^01 + \,k (fffcl^ll + 5fe2^2l)] 



5fc2-''^23 

-1 

(3.8c) 



Our notation applies to a lattice made of N sites, peri- 
odic boundary conditions are assumed, /? is the inverse 
temperature (the Boltzmann constant is set to unity), 
finally q = (itu^ , q) and fc = (iw„,fe) are three vectors 
with bosonic and fermionic Matsubara frequencies, re- 
spectively, while q and k belong to the first BZ. It is 
straightforward to modify this notation for the case of 
the continuum limit. 

After performing the summation over the fermionic 
Matsubara frequencies in Eq. ( |3.8a ) and with some addi- 
tional lengthy algebra, the dynamical Pauli susceptibility 
tensor simplifies to 



1 



yXau) 



AN ^ H i^M^ 

fc A, A', 7, 7' 



A, A' 



q \ / fc,q 

UT>{lEux)-UT>iiE^+^w) 



k+qX' 



IVJ, 



(3.9a) 



where 7,7' = ± and n^v = 0,1,2,3, the single-particle 
dispersion in the superconducting state E^,-^ is defined in 
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Eq. (2.13) for the helicities A = ±, while 

1 



/fd(^) 



1 



(3.9b) 



is the Fermi-Dirac function. The vertex is given by [g^^ 
9kl\9k\) 



00 
60 

Od 

^A,A' 
bd 



k.q 



k.q 



k.q 



k^q 



l + ^^'Sfc Sk+q: 

- ^' Sb-k+q + ^9b-k + ^^'i£ahcffa;fc5c;fc+qi 
^' 9d:k+q + ^9d:k ^ '^'^adcfj a-k9 c;k+qi 

- ^bd + 9a-kfbd9c-k+q 

+ i^abd {^9a;k ~ ^' 9a:,k+q) ' 



(3.9c) 

where the tensor /^^ is defined by 

fbd -5,,^,, - 5,Ad + ^ad^bc - Pdb- (3.9d) 
FinaUy, the coherence factors are given by 

i£,k,x£,k+q,\' — ^k.X^k+q,\' 



' ^A,A',7,7' 
. "-"00 



k.q 



= 1 + 77 



Ek+q,\'Ek^\ 



A, A ,7,7' 
60 



_ / ^A,A',7,7'\ _ Cfc,A / 'Cfc+q.A' 



A. A' .7,7' 
hd 



k,q 



■77 



-Bfc^A Ek+q,\' 
i£,k,x£,k+q,\' + ^k,\^k+q.\' 



Ek+q.y Ek.X 



(3.9e) 



where the single-particle dispersion in the normal state 
)^ and the superconducting pair potentials are 
defined in Eq. ( ^.4bD and in Eq. ( ^.9a ), respectively. 
The dynamical Pauli susceptibility of the normal state 
is obtained by taking the limit A^ a ^ 0, A = ±, 
supplemented by the replacements ^Ej^^ Ep,^ and 
7'-E'fc_l_qA' ~^ -^fc+qA' to remove the particle-hole 

symmetry. 

We then recover the result from Ref. bll If we fur- 



thermore set = in Eq. (3.9a), we obtain the static 



susceptibility of the normal state. At the Rashba-Dirac 
point, i.e., for £,1^^ — X\gi,\, the following components of 
the static susceptibility vanish: Xoi> X02, Xio, X20' X13, 
X23> and X32. 



We now consider the isotropic continuum limit (2.1) 
(the Rashba-Dirac limit) together with an attractive con- 
tact density-density interaction — |y|(S(r — r'), which in- 
duces purely singlet pairing Ag as we will show in Sec. 
IV A. The gap equation at the chemical potential yU and 



the inverse temperature /3 is 



deu{e) tanh/3£'^(e)/2 



(3.10a) 




FIG. 2: (Color online) Temperature dependence of the 
selfconsistent superconducting gap (red, normalized by the 
value at zero temperature) and the in-plane and out-of-plane 
Pauli magnetic susceptibility (red, normalized by the max- 
imum of the out-of-plane susceptibility) for ^ = Scjjj and 
jV|!/(ajD) = 1. Here, cjd is the Debye cutoff used for the gap 
equations. 



where the "Debye" energy cutoff oj^ has been introduced, 
E^{e) := ^(F^^vTTAf, (3.10b) 

and 



v{e) 



(3.10c) 



is the Rashba-Dirac density of states per unit area. The 
temperature dependence of the static Pauli magnetic sus- 
ceptibility for an out-of-plane uniform applied magnetic 
field is then given by X33 = 2xii = 2^22 with 



X33 « XI 

A=± 



A^(r) 



(27r)2 
1 



^fcA 



?fc,A Sfc, 



■ tanh 



(3.11) 



fcA 



We plot the temperature dependence of the self- 
consistent pair potential Ag(r) and of Xz^iT) in Fig. |^. 
First, we note that XzaiT = 0) 7^ 0. |-Tl|iis is a direct 
consequence of the spin-orbit coupling.llII^E2l Second, we 
note that X33 {T) decreases as a function of temperature 
beyond the critical temperature, i.e., when T > T^. Al- 
though the finite value of X33 (2^r-^ Q^nis typical of noncen- 
trosymmetric superconductors JiH^El X33(^ > ^c) only 
saturates to a value proportion al to the density of stat e 
at the Fermi level in the regime ( |l . l|) . In the regime (^]^) , 
the decreas e of X33 (T > Tj, ) can be understood with the 
help of Eq. ( |3.5[) if ^ is sub stitu ted for T to mimic thermal 
population. Indeed, Eq. ( |3.5| ) implies that the normal- 
state X33(m) at T = increases with fjL if the Fermi level 
is below the Rashba-Dirac point, but decreases with ^ if 
the Fermi level is above the Rashba-Dirac point, for the 
states above the Rashba-Dirac point give a contribution 
that cancels part of the susceptibility coming from the 
states below the Rashba-Dirac point. 
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IV. SUPERCONDUCTING INSTABILITIES 



Density-density interaction 



To obtain the BdG Hamiltonian ( 2.5a ) self- 
consistently, we consider first a density-density in- 
teraction given by 



q fc,s=t.-l 

where is an even function of momentum. Normal 
ordering yields 



^ E ^'^k+qs'^ks (4.1) 



ks 



kM' s,s' 



(4.2) 



After renormalization of the chemical potential and re- 
striction of the normal-ordered interaction to the scat- 
tering of Cooper pairs with vanishing centcr-of-mass mo- 
mentum, wc obtain the reduced Hamiltonian 



^v'^ — 9 E ^k-p E ^^ks^^-ks'^-ps'^^ps 



(4.3) 



kyP 



We show in Appendix ^ that the reduced interac- 
tion (|4.3|) has the representation 



H 



red 



JE E yk-pe<-^-^^) 

k,p A.A'=±1 

X [cos {(fip - ifk) + W] aLa-fcAa-pA'OpA' 



(4.4) 

The terms • • • that were omitted involve pairs of creation 
or of annihilation operators of opposite helicities. We 
ignore these terms because we are going to perform a 
mean- field approximation for Cooper pairs made of tim e- 
reversed helicity single-particle states from Eq. ( 2.7c ). 

We define the mean-field superconducting order pa- 
rameters to be 



-•-fcA 



(4.5a) 



where A = ±. The angular bracket represents the sta- 
tistical averaging at inverse temperature j3 and chemical 
potential ^. We also define the mean-field helicity pairing 
potentials to be 



^kx — J E ^k-p \-^^' (^P " "fk) + 1] Ky 



p,\' 



(4.5b) 



fcA 



where A = ±.Thc mean-field superconducting order pa- 



rameter (4.5a) and the pairing potentials (4.5b) enter the 



(mean- field) BdG Hamiltonian of the form ( 2.12| ) and 
obey the self-consistent conditions 



"pA 



2E. 



tanh(/3£;pj2) 



pA 



(4.6a) 



where the single-particle dispers ion i n the superconduct- 
ing state is defined in Eq. ( ^.13| ). 

If the pairing interaction is independent of momentum 
(i.e., a contact interaction in space ), the summation over 
p on the right-hand side of Eq. ( 4.5b ) cancels the de- 
pendence on A. Hence, both order parameters are then 
equal A^.^ = and wc can sec from the transforma- 
tion (2.10b) that the pairing potential will be of pure 
spin-singlet nature. Observe that this result is indepen- 
dent of the noninteracting p art of the Hamiltonian, and 
thus valid for both models (2.1) and (2.2). It was also 
found in the context of 3D noncentrosymmetric super- 
conductors in Ref. 

We have also solved self-consistently the gap equation 
with the Dirac dispersion (e^ = 0) for a pairing inter- 
action that is isotropic in momentum space = V^^y 
When the chemical potential is much larger than the 
transition temperature 3> 1, we have found that 
the triplet component never exceeds the singlet compo- 
nent of the superconducting pairing potential if the pair- 
ing interaction V|g| never changes sign as a function of 
\q\. The latter is true for most of the commonly used 
model interactions, except Cooper pairing mediated by 
the Friedcl oscillations induced by the screeaiag of the 
Coulomb repulsive interaction, for example. ljEZI 

The density-density interaction as considered here 
might provide a model for the pairing interaction re- 
cently discovered at the superconducting interfaces in 
LaAlOa/SrTiOg (Ref. ||) or in some electrolyte/SrTiOa 
(Ref. ^9|) that feature a low density and high mobility of 
the charge carriers. 



B. Heisenberg interaction 

As a second example, we study the SU(2) preserving 
spin-density-spin-density interaction 



2 E ■^q'^q ' 2 E ^k 

q k\s,s 



k+qs'^s.s'Cks 



(4.7) 



where is an even functio n of m omentum. Proceeding 
in the same way as in Sec. [VA, we obtain the reduced 



Hamiltonian for the scattering of Cooper pairs with van- 
ishing ccnter-of-mass momentum 



- ^ E E ^fc- 



P" Sl^«4 



(4.8) 
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As is shown in Appendix ^ the reduced interaction (4.8 
has the following representation in the helicity basis 



k,p A,A' = ±1 



X [cos {(pp - - 3AA'] aL«-fcA«-pA'apA' 



(4.9) 



The terms • • • that were omitted, just as in Sec. |va| 
where we studied density-density interactions, involve 
pairs of creation or of annihilation operators of opposite 
helicities. We ignore these terms because we are going to 
perform a mean-field approximation with pairs of time- 
reversed helicity single-particle states from Eq. (2.7c). 

Again, we define the mean-field superconducting order 
parameters to be 



^feA := '^e"^''° {a--kxa-kx)p,^ = 



-fcA 



(4.10a) 



where A = ± and the angular bracket represents the sta- 
tistical averaging at inverse temperature (3 and chemical 
potential fj,. The mean-field helicity pairing potentials is 
defined to be 



p,A' 



(4.10b) 



fcA 



where A = ±. 

Together with the superconducting order parame- 
ters (4.10a), they obey the self-consistent conditions 



( |4.6aD . 

The term proportional to A represents the triplet com- 
ponent of the pairing potential while the constant term 
in the square bracket gives the singlet component. 

We have solved self-consistently the gap equation with 
the Dirac dispersion (i.e., = 0) for a pairing interac- 
tion that is isotropic in momentum space Jq ~ J\q\- 
with the case of the density-density interaction, we have 
found under the assumption |^/3| ^ 1 that the triplet 
component never exceeds the singlet component of the 
superconducting pairing potential if the pairing interac- 
tion J|g| never changes sign as a function of \q\. 

The Heisenbcrg interaction is attractive (repulsive) in 
the spin-singlet channel and repulsive (attractive) in the 
spin-triplet channel for > (J^ < 0). In centrosym- 
metric superconductors, this property leads the way to- 
ward spin-fluctuation mediated spin-triplet superconduc- 
tivity for Jq < 0. If inversion symmetry is broken, how- 
ever, spin-singlet and spin-triplet pairing channels are not 
separated and for the case Jg < 0, the interaction is alto- 
gether not attractive. Hence, the Heisenberg interaction 
will not lead to triplet (dominated) superconductivity in 
the same fashion as in centrosymmetric superconductors. 



C. Dzyaloshinskii-Moriya interaction 

Finally, we study a spin-density-spin-density interac- 
tion of Dzyaloshinskii-Moriya type, which requires the 
breaking of inversion symmetry to be present. Let the 
coefficient be a three vector with vanishing z com- 
ponent for our case. It shares the symmetry of g^, in 
particular, it is odd under q — > —q. The Dzyaloshinskii- 
Moriya interaction interaction is then 



H 



DM 



(4.11) 



The reduced Hamiltonian for the scattering of Cooper 
pairs with vanishing center of mass momentum reads 



rrred 
-"DM 



^E E ^fc-p-K, 

^ '^ks,'^-ks„^-ps2'^ps^- 



A a. 



(4.12) 



As is shown in Appendix ^ the reduced interac- 
tion (4.12) has the representation 



rrrod 
^DM 



^E E ^'^^-^^^ 

k,p A,A'=±1 

^ p ■ (^•^j^ ^ ^ "feA«-feA"-pA'apA' 



(4.13) 

Once again, the terms • • • that were omitted involve a 
pair of creation or of annihilation operators of opposite 
helicities while we keep only pairs mad e of t ime-reversed 
helicity single-particle states from Eq. (2.7c). 

The mean-field superconducting order parameters are 
again defined to be 



kX^kX/ p^fj 



kX 



(4.14a) 



where A = ± and the angular bracket represents the sta- 
tistical averaging at inverse temperature /3 and chemical 
potential fi. We also define the mean-field helicity pairing 
potentials to be 



A 



fcA 



—ol^D,_p (^A — -A— 



p,A' 



-'pA' 



(4.14b) 



= A 



-fcA 



where A = ± . Tog ether with the superconducting order 
para meters (4.14a), they obey the self-consistent condi- 
tions ( 4.6a ). 

We have solved self-consistently the gap equation with 
the Dirac-dispersion (e^ = 0) for the model interaction 
I , where a and D are parameters. 



We have found that, depending on the chemical poten- 
tial and the cutoff parameter a, the triplet component 
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can exceed the singlet component of the superconduct- 
ing pairing potential. This is in contrast to the results 
from the density-density interaction and the Heisenberg 
interaction, where the singlet component is dominant. 
However, as the Dzyaloshinskii-Moriya interaction arises 
in second-order perturbation theory from an exchange 
interaction, it should not be considered on its own. 



D. 



Superconductivity with in-plane magnetic field 



assume that only electrons of a single helicity A form 
Cooper pairs. From here on, we denote the center-of- 
mass momentum of Cooper pairs by q while k and k' 
refer to the relative coordinate of the paired electrons. 
We assume that a single wave vector q for the modula- 
tion of the pairing potential will be energetically favor- 
able, rather than a set of degenerate wave vectors. With 
these simplifications, the self-consistent gap equation for 
the pair potential Kkx[q) at temperature T close to the 
superconducting transition temperature and for TV sites 



A well established result for 2D noncentrosymmet- 
ric su perco nductors with the Rashba spin-orbit cou- 
pling (2.2c) is that the superconducting pair potential 
acquires a real-space modulation in the pre senc e^ a Zee- 
man coupling to an in-plane magnetic field.llZl&Ea An in- 
plane magnetic field shifts the Fermi sea away from the 
center of the Brillouin zone and as a result, the electrons 
with opposite wave vectors are not degenerate in energy 
anymore. 

A similar effect is expected for the 2D Rashba-Dirac 
model subject to this study. We shall demonstrate this 
for a momentum- indepe ndent density-density interaction 
as was discussed in Sec. LVA. 



The noninteracting Hamiltonian (2.4a) in the Rashba- 
Dirac limit, i.e., for e^. = 0, is altered in the presence of 
an in-plane magnetic field 



according to 



(4.15) 



n 



0;fc 



''■O-.k 



keBZ keBZ 

B) cr, 



-Mc^o + idk 

e 



(4.16) 



The single-particle dispersion is now given by 
i^^ = -^^±\g^-B\. 



(4.17) 



Accordingly, the phase factor entering the transforma- 
tion ( 2.4c ) between the laboratory basis and the helicity 
basis is changed to 



3ki 



iSo 



A pair of electrons on the Fermi surface without mag- 
netic field with opposite wave vectors k and —A; is not 
degenerate in energy anymore in the presence of B, for 



'^^k'+q/2 '^^k+q/2 



k' 



X COS (^l±lll_f?±^\ A^,^^^^ 



X G 



(0) 



(0) 



fc'+q/2,iLJ^;A —k'-\-q/2, — iuj^^;X 



(4.20) 



where the single-particle Green's function in the normal 
state for electrons with helicity A = ± is given by 



1 



(4.21) 



For s-wave pairing, the pairing p otent ial Ak\{q) is inde- 
pendent of fe. The gap equation ( 4.21 ) simplifies to, after 
performing the summation over Matsubara frequencies. 



1 = cos 



'P^k'+q/2 V^k+q/2 



xcosp:±^zV^ /-,^,. 



(4.22a) 



with the function 
„ tanh ■ 

fB — 
J k.q;X ' 



2T 



tanh ■ 



2T 



fc+(q/2)A 



^k+{q/2)X 



(4.22b) 



In the Rashba-Dirac continuum limit (2.1), the disper- 
sion (4.17) together with the transformation (4. IS) and 
the thermal factor (4.22b) obey the symmetries 



cB 

^k+(qj2)\ 
B 

'^k+{qj2) 



pB= 
^k\ 



B=0 
fcA I 



(4.23) 



/i 



fc,(qo/2);A 



fc.O;A' 



k± 



cB 
?-fe± 



T2 



9k - B 

\9k\ 



O 



\9k\ 



(4.19) 



It might thus be energetically more favorable to pair 
electrons with the same energy but with a finite center- 
of-mass momentum, than to pair electrons with vanish- 
ing center-of-mass momentum. For simplicity, we also 



with Qg = 2B A e^Uhv-^-o) being proportional to the 
shift of the Fermi surface induced by B. Hence, the 
gap equation for the superconducting condensate with 
the center-of-mass momentum qg in the presence of the 
in-plane magnetic field B is the same as the gap equa- 
tion in the absence of any in-plane magnetic field for 
a condensate with vanishing center-of-mass momentum. 
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A condensate with vanishing center-of-mass momentum 
has the largest transition temperature. Hence, we deduce 
from the symmetry (4.23) that a superconducting order 
parameter with the center of mass momentum Qq nucle- 
ates in the presence of an in-plane magnetic field. The 
wave vector of the modulated pairing potential is per- 
pendicular to the magnetic field in the plane and is in- 
dependent of the chemical potential. It also follows that 
the critical temperature of superconductivity is not sup- 
pressed by the magnetic field in the Rashba-Dirac con- 
tinuum limit ( p.ip and by the consideration of only one 
helicity. 

In contrast to this simple result, the center-of-mass mo- 
mentum of Cooper pairs in 2D noncentrosymmetr ic su - 
perconductors with the Rashba spin-orbit coupling ( 2.2c ) 
is selected by a delicate energetical compromise on how 
the two helicity-resolved Fermi surfaces are shifted in op- 
posite directions in momentum space. 



V. MEAN-FIELD PHASE DIAGRAM 

It is time to explore the mean- field phase diagram that 
follows from Hamiltonian ( ^.lOa ) in the parameter space 
spanned by the choice made for the normal-state disper- 
sion and for the pair potentials. To this end, we consider 



the parameter space spanned by the couplings fi, A^, 
and entering the mean-field Hamiltonian. A mean- 
field phase corresponds to a connected region in parame- 
ter space characterized by a nonvanishing gap. We shall 
then introduc e in Se c. |Vl| point defects in the mean-field 
Hamiltonian ( 2.10a ), i.e., superconducting vortices with 
unit circulation, and compute the parity of the number 
of zero modes they bind to characterize the topological 
nature of the mean-field phases separated in parameter 
space by gap-closing boundaries. 

For convenience, we reca ll that the BdG Hamiltonian 
is, in the $ representation (2.10a), 



A. 



.k 



t.k' 

A„ 



A 



S.fc 



A. 



A 



A 



S.fc 



where the normal-state dispersion is specified by 



^k = ^-k £ M £ IK, 

Ak^\9k\ = ^-k>0, 

= arctan^^ e [0, 2n[, 

9k;l 



9k = 9 k e 



(5.1a) 



(5.1b) 



the singlet-pair potential A^ and the triplet-pair po- 
tential Aj. transform according to any trivial irre- 
ducible representation of the space group consistent with 
A^. fcC^"'''" being single valued. In the isotropic contin- 
uum limit, we thus assume that the singlet-pair poten- 
tial is constant while the triplet-pair potential A^(fc) can 
be factorizcd into a real number A^ times some strictly 



increasing positive function / with at least a first-order 
zero at the origin such that (i) it saturates to unity for 
large positive argument and (ii) is invertible on the pos- 
itive real axis with the inverse /~^, say, for instance, 
f{x) := tanhx, i.e.. 



A,(fc) = A,/(|fe|/fcJ 



(5.1c) 



for some wavelength k^. > that defines the size of the 
core of the vortex exp [ — i<y3(fc)] at the origin in fc space. 
The aim of this section is to identify when the quasipar- 
ticle spectral gap vanishes as a function of the parame- 
ters Ag , A^ , and /i for a given dispersion relation in the 
isotropic continuum limit, i.e., we need the eigenvalues 
of Eq. ( |5l| ). 

To this end, we first square both sides of Eq. ( ^.l| ), 
finding the block diagonal form 



fc 

A, 



2 [(^fc - m) ^fc + W,k] cos ifik 
2 [(^fc - Ai) ^fc + W,k] sin fk ^2- 



(5.2) 



The four eigenvalues of V.^ are 



^fc;A,± =±\/(£fc-/^ + ^^fc) 



AA 



t,fc / ' 



(5.3) 



where A = ±. All nonvanishing energy eigenvalues come 
in pairs with opposite signs. This spectral symmetry is 
a consequence of the particle-hole transformation [X^^^ 
was defined in Eq. 1 0^ )] 



^22T~L-k^22 — ^T~Lk- 



(5.4) 



The Hamiltonian "H^ also features a helical symmetry 
given by 

(5fcl-''^10 + 9k2^2Q) T~Lk {9kl^lQ + 9k2^2o) = '^fc- (5-5) 

Viewing the Rashba spin-orbit coupling as a fictitious 
magnetic field along a fc-dcpendcnt direction, the helical 
symmetry (5.5) reflects the conservation of spin along 
this direction in momentum space. 
For completeness, TRS is nothing but 



^2o'^-k^20 



(5.6) 



in the basis of Eq. (2.51:) that wc have chosen. 

Zero modes are vanishing eigenvalues of T-Lk, i-e., they 
arc the solutions to 



= det Uf, 



(£fc -M + ^fc) 



A 



t.fcy 



(5.7) 



(£fc - M - Af,f + (A, 



H,k 
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FIG. 3: (Color on line) Mean- field phase boundary in the 
Rashba-Dirac limit (IS.ltj). 



There arc thus two possibilities to get zero modes. Either 

case(+): = + = A. + Aj^^, (5.8a) 

or 

case(-): = £fc-/i-^fc, = A^ - A^ (5.8b) 



Equation (|5.8a| ) requires that the A = + heheity gap 
vanishes on the A = + helicity Fermi surface. Equa- 



tion ( 5.8b ) requires that the A = — helicity gap vanishes 
on the A = — helicity Fermi surface. The condition 



M + XAf, = 



(5.9a) 



on the normal-state dispersion determines the Fermi sur- 
faces 

FS^:={k\e^~fi + XA^=Q}. (5.9b) 
The condition 

A, + XA, 4 = (5.10a) 

on the pairing potentials determines the momenta for 
which the superconducting single-particle gap vanishes 

SC^:={fc|A,±At4=0}. (5.10b) 

Conditions {+) or ( — ) are satisfied along the points 

FS;, n SC;, 7^ 0, A = ±. (5.11) 

(in other words, the Fermi surfaces cross the supercon- 
ducting single-particle nodal surfaces). 

A. Isotropic continuum limit 

We work in the continuum limit with the upper bound 
A and the lower bound —A to the single-particle mean- 
field spectrum, as is appropriate for the surface states of 
a 3D band insulator. We assume that the SRS dispersion 



e, the Rashba dispersion A, and profile / of the vortex 
exp [ — iv5(fc)] at the origin in k space are smooth func- 
tions of |fc|. For the analysis to come, it is useful to define 
the dimensionless quantity 



k^|fc|/fc, 



(5.12) 



We define the 2D parameter space with A^/A^ as the 
horizontal axis and fi as the vertical axis. For any finite 
positive singlet pairing potential Ag ^ 0, we show 

1. That there are two nonintersecting curves (to sim- 
plify the notation A — > oo) 



Ai+ : (-c»,-l] ^M, 

A,/A, ^ A'+(At/AJ, 

Ai_ ; [l,oo) — > R, 

A,/A3^A^_(A,/AJ, 



(5.13) 



defined by the condition ( ^.8a ) for A = + _and ( ^.8b| ) 
for A = — at which the BdG spectrum ( ^.3[ ) is gap- 
less. 

2. The curves fi^ are one-to-one reparameterizations 
of the dispersions C;^(|A;|) with A = ±. 

3. How the two curves fi^ change upon changing the 
topology of the Fermi surfaces. 

For the superconducting single-particle gap to vanish, 
we must choose 



As 

A = -sgn^ 



(5.14a) 



in Eq. ( ^.lOa ) from which the implicit definition 

A. 



< /(k) 

of k follows. Hence, k is the function 
k : (l,oo) — > R, 



(5.14b) 



At 




\ 




A. 


) 


A. 




A. 


At 





(5.15) 



which is not defined whenever the superconducting 
single-particle gap does not close, i.e., when \A^\ < |Ag|. 
Claims |^ and || follow with the definition 



/iA(k) :=£(k) + AA(k) 



(5.16) 



where k and A were themselves defined in Eq. (5.14) and 
with the momentum core size taken to be unity. 

To illustrate how the topology of the normal-state dis- 
persion changes the curves ij,^ with A = ±, we make the 
(electronlike) parabolic approximation 



e(fc) 



h^\k\- 

2m 



m > 0, A{k) = hvj^j^\k\, (5.17) 
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and choose the momentum-vortex profile 

/(a;)=tanhx (5.18) 

with the momentum core size taken to be unity. We 
consider the Rashba-Dirac hmit 

m = oo (5.19) 

first, as is ihustratcd in Fig. ^. The two curves 
/x_l_(zV/Ag) where the gap vanishes are indicated in 
Fig. H. They arc obtained via the rcparameterization 
of the dispersions ^_|_(|A;|), as depicted on the insets on 
the second and fourth quadrants. There is a one-to-one 
correspondence between the thick fines in the insets and 
the curves /x_(_(At/Ag). 

In Fig. ^, we look at the mean-field phases that can be 
identified given the gap-closing curves /^_|_(A^/Ag). Here, 
one must notice that taking A^/Ag — > oo, for a given 
chemical potential such that the gap does not close, is 
connected to the path originating from A^./A^ — >■ — oo. 
For instance, one can send Ag — >• so that it changes sign 
while holding Aj constant but with a given /i such that 
the gap does not close. Therefore, the regions depicted in 
Fig. raare connected upon folding the horizontal axis into 
a circle (the plane into a cylinder) . Then, because of the 
topology of the curves /x_(_(Aj/Ag), any one region can 
be connected to any other without crossing these curves, 
and hence there is a single phase for the system, which 
we denote by I. 

For any finite curvature of the dispersion e(fc), i.e.. 



< m < oo. 



(5.20) 



we find the boundaries shown in Fig. |[ Again, the two 
curves fi^{A^/ A^) where the gap vanishes are obtained 
via the rcparameterization of the dispersions ^_|_(|fc|), as 
depicted on the insets on the second and fourth quad- 
rants. We see that the topology of the curve /x_|_(|At/Ag|) 
that tracks the normal-state dispersion ^_|_(|A;|) is insen- 
sitive to tuning m from infinity to any finite value. This 
is not so for the topology of the curve /^_(|Aj/Ag|) that 
tracks the normal-state dispersion ^_(|A;|). This curve 
is dramatically infiuenced by the nonmonotonous depen- 
dence of ^_(|fc|) on |fc| for any finite curvature, i.e., any 
mass m < oo. In the Rashba-Dirac limit m = oo, 
/i_(|Aj/Ag|) is strictly negative, and ^_(1^) — oo. 
But when m is finite, /x_(l_|_) — >■ oo. 

The distinct phases in the m finite case are identified in 
Fig. § If the regions with A^/A^ ±oo are identified, 
then the two regions II and III are always separated from 
each other by the two curves fj,^{\A^/ AJ). There is no 
path connecting the two regions II and III without ever 
closing the mean-field spectral gap: this is a necessary 
(but not sufficient) condition for these to be two distinct 
phases. 
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FIG. 4: (Color online) Mean-field phas e bou ndary away from 
the Rashba-Dirac limit, i.e., when Eq. (5.20) holds. 



B. Anisotropic case 

The boundaries in the A^-fj. plane at which the BdG 
single-particle spectrum closes for an anisotropic con- 
tinuum dispersion or for a two-dimensional lattice are 
more difficult to determine. Indeed, a technical difficulty 
brought about by the loss of continuous rotational sym- 
metry is that it is not possible anymore to characterize 
the nodes of the normal-state dispersion or the nodes 
of the superconducting gaps with a single wave number. 
This could result in those boundaries acquiring a thick- 
ness (i.e., a finite area).E2l 

For a 2D lattice model, a qualitative difference with the 
continuum limit that is of no relevance to this section is 
the fermion doubling and its consequences for the exis- 
tence and the stability of Majorana fermions in the core 
of superconducting vortices. This is the subject of the 
ensuing section in which we search for Majorana modes 
in the core of defects (vortices) of the superconducting 
order parameter and we probe their stability under adia- 
batic changes of the bulk parameters (i.e., far away from 
the vortices). 



VI. MAJORANA FERMIONS 

Caroli et al. showed in Ref. ^ that isolated vor- 
tices in a weakly coupled type II s-wave superconductor 
with TRS and SRS support a discrete set of finite-energy 
bound states with a level spacing of order of the ratio 
of the squared single-particle bulk superconducting gap 
to the bandwidth. There is no bound state at the Fermi 
energy bound to the core of vortices in this case. 

Jackiw and Rossi showed in Ref. ^ that, in two 
space and one time dimensions quantum electrodynamics 
(QED2+1) coupled to one scalar Higgs field, an isolated 
static defect in the Higgs field, i.e., a single vortex with 
vorticity N, supports N bound states that arc all pinned 
to the zero energy. These N bound states arc N Ma- 
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jorana fcrmions. Aif-index theorem for this result was 
proved by Weinbeijg£3 

Read and GreencHl showed that a two-dimensional chi- 
ral ± ipy superconductor supports a Majorana mode 
bound to the core of an isolated half vortex. 

We are going to show that (i) an isolated vortex with 
unit vorticity in the singlet-pair potential binds a single 
Majorana mode in region I of Fig. |[ (ii) an isolated vor- 



tex with unit vorticity in the triplet-pair potential binds 
two Majorana modes in region II of Fig. and (iii) iso- 
lated vortices in region III of Fig. || do not bind Majorana 
fermions. We will start by reviewing the derivation of the 
Jackiw-Rossi Majorana mode that applies to region I of 
Fig. ||. We will then discuss region II and III in Fig. ^. 

We work in the isotropic continu um lim it with the 
Hamiltonian in the $ representation ( [2.10a ) given by 



hv^^k sik,-k) i{A,(z,z),/_(M)^| 

H.c. H.c. ~e{k,k) 
y H.c. H.c. —hv^j^k —e{k,k) j 



-hv^^k 
— e(fc, fc) 



(6.1a) 



Here, the SRS normal-state dispersion is parabolic 



in terms of which 



£(fc, fc) 



h^\k\' 

2m 



(6.1b) 



where the real valued /i is the chemical potential. More- 
over, the singlet-pair potential Ag(z, z) has a unit vortex 
at the origin of the complex-z plane with the character- 
istic core size and saturates to the bulk value for 
|z| » 4: does the triplet-pair potential Aj(z, z) with 
the characteristic core size and the bulk value A^ for 
|z| ^ l^. The bulk values Ag and A^ share a common 
phase that can be removed by a global gauge transfor- 
mation up to a relative sign. The function / that guar- 
antees single valuedness of the Hamiltonian was defined 
in Eq. (5.1c). The anticommutators in the antidiagonal 
are needed since translation invariance has been broken. 
We are using the complex notation 



k . — Aj-^ ^^ ^^2' ^ * — ^1 ^^2' 



z :— 1Z2, 



together with the algebra 



(6.2a) 



(6.2b) 



K>y = i<5a&: a,&=l,2, 
or, equivalently, 

[z, fc] = [z, fc] = 2i, [z, fc] = [z, fc] = 0. (6.2c) 



A representation of the algebra ( 6.2c ) is given by 

k = -219^, k = -2i9^. (6.2d) 
The representation dual to Eq. ( |6.2d ) is 

z = 2i9^, z 219^. (6.2e) 
We shall rely on the polar coordinate representation 

k = Ke+'^ (6.3) 



z = 2ia^. = ie+'^ ( 5, + -d^ ) , 



z = 219;. = ie 



9. - -a.. 



(6.4) 



We choose to represent z, z as differential operators of 
functions of fc, k instead of the other way around because 
this can bring a simplification in the solution of the zero 
modes. By solving for the wave functions of the zero 
modes in momentum space, we take advantage of the fact 
that we have a first-order differential equation instead 
of a second-order one, which would be the case had we 
chosen to solve for the wave functions in position space. 
This simplification works because we deform the profile 
of the vortex without changing the fact that the solutions 
are precisely at energy E' = 0, as we discuss below. 

Instead of facing the difficulty to solv e an alytically for 
the spectrum of the BdG Hamiltonian (S.l), we are thus 
going to make approximations that are motivated by the 
mean-field phase diagram of Sec. |^. 

We are first going to take the Rashba-Dirac limit at 
/i = (the Rashba-Dirac point) without triplet-pair po- 
tential, Aj = 0. This is nothing but the origin of region 
I in Fig. [3. In this limit, Hamiltonian (6.1) is the direct 
sum of two 2x2 Hamiltonians. 



We are then going to take the Fermi limit 







with fi > without singlet-pair potential, A^ = 0, i.e., 
the vertical half line at i nfin ity in region I of Fig. m In 
this limit, Hamiltonian (3.1) is again the direct sum of 
two 2x2 Hamiltonians. 

Even after these simplifications, the spectrum with a 
vortex at the origin is difficult to compute. One must 
solve a system of two coupled partial differential equa- 
tions that depends on the non universal details encoded 
by the profile of the vortices in the superconducting pair 
potentials and by the profile / for the vortex in the d 
vector. This microscopic information does influence the 
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FIG. 5: (Color online) Qualitative comparison of the bound 
state spectra of a deep and narrow vortex (left) and a wide 
and shallow vortex (right). While the former supports only 
few finite-energy bound states (blue), in the spectrum of the 
latter more states are allowed. However, the existence of a 
zero-energy mode (red) is independent of the details of the 
regime. 



discrete spectrum that represents the bound states with 
nonvanishing energies attached to the core of the vortex. 
For example, if the profile of the vortex is deep due to 
a large bulk gap and narrow due to a small character- 
istic size, there should be very few bound states with 
nonvanishing energies (see Fig. On the other hand, 
in the opposite limit of a shallow and smooth profile for 
the vortex, many bound states with nonvanishing ener- 
gies extending far away from the vortex core are to be 
expected (see Fig. ||). 

However, we arc not after the full spectrum of states 
bound to the core of a vortex. We are only seeking 
the conditions under which Majorana bound states, i.e., 
bound states pinned at the normal-state chemical poten- 
tial, are present. The very existence of a Majorana state 
does not depend on the profiles of the vortices in real 
and momentum space as long as the nonvanishing energy 
spectrum of bound states remains discrete and separated 
from the zero energy. This suggests choosing the vortex 
profile 



3.5a) 



for the singlet (s) or triplet (t) component of the pair 
potential in real space, respectively, and the vortex profile 



3.5b) 



for the d vector in momentum space. This approxima- 
tion has the merit of linearizing the spectral eigenvalue 
problem. 



In this limit and when the singlet-pair potential vanishes, 
the bulk gap closes so that vortices in the triplet-pair 
potential are ill defi ned. 

In the limit ( |6.6| ) and when the triplet-pair potential 
vanishes, after setting ?i = Wpj-, = 1, 



/ k 
k 
H.c. 

\ H.c. 



\ 





-k 
/ 



(6.7a) 



decomposes into the direct sum of the 2x2 Hcrmitian 
operators 



and 



^(2) 



-k 



' _k 

-iT-Z 



(6.7b) 



(6.7c) 



Jackiw and Rossi showed that Hamiltonian (6.1a) 
with a vortex in Ag(z,z) satisfying < oo and 
I Ag(z, z)| i^i^oo < CO- supports one and only one bound 
state and that this bound state is pinned to the chemical 
potential, i.e., to the midgap of the BdG Hamiltonian. 
We will show explicitly that Hamiltonian (6/7), with an 
unbounded vortex profile, also has a singly degenerate 
solution, and thus there is one and only one Majorana 
fermion in the Rashba-Dirac limit at the Rashba-Dirac 
point if the pair potential is pure singlet. The stability 
of this Majorana fermion away from the Rashba-Dirac 
point or in the presence of a triplet-pair potential is guar- 
anteed by the fact that the particle-hole symmetry of the 
eigenvalue spectrum precludes the migration of the zero 
mode as long as the gap does not close. For large val- 
ues of the chemical potential, it is natural to anticipate 
that many more bound states will have peeled off from 
the continuum with a level spacing a la Caroli-dc-Gennes. 
This expectation is consistent with the computation from 
Rcf. |6^ of the states bound to a vortex as a function of 
the chemical potential for pure-singlet superconducting 
graphene. 



Proof. Let 



(6.8) 



A. The Rashba-Dirac limit 

The Rashba-Dirac limit is defined by the condition 

m = oo. (6.6a) 

The Rashba-Dirac-point limit is defined by the additional 
condition 



^ = 0. 



3.6b) 



We seek the solutions to 

= fcu*^) +czv'^^\ 
= czu^^^ - kv^^\ 



and 



= fcw^^) +czv^'^\ 







,(2) 



kv^ 



(6.9a) 



(6.9b) 
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respectively. If we take the c omple x con jugate of the 
second condition in Eqs. ( 3. 9a ) and ( 6.9b ), respectively, 
we get 



and 



= fc?i(i' + czw(i), 
= -fcw(i) +czu(i' 

= h^^'> + czv^'^\ 
= -fci;(2' + czvP-'^ 

respectively. For any j = 1,2, if 



is a zero mode, so is 



± 



^u(j)(-fc,-fc) 
Hence, we try the ansatz 



$2^(/c,/c) := 



u 



U) 



{k~k) 



^±u(j)(-fc,-fc), 
where 

= ku'^^\k, k) + (±)czu(i)(-/fc, -k) 
and 

fcu(2)(fc, k) + (±)czw(2)(_fc^ 



We choose a gauge in which 

c := ic 

is real and make the ansatz 



where the real- valued (k) satisfy 



K + ( 9„ + - 



and 



0^[n+{±)~cd,]w'i\K), 



(6.10a) 



(6.10b) 



(6.11) 



(6.12) 



(6.13a) 



(6.13b) 



(6.13c) 



(6.14) 



(6.15a) 



(6.15b) 



(6.15c) 



respec tively. The formal solutions to Eqs. (6.15b) and 
( |6.15c ) are 



w^^' (k) = w^^' (kq) exp < 



(1) 



dn' 



C K 



(6.16a) 



and 



(2), 



respectively. Only 



hi 



,.2 ^2 
^0 



sgncv"/ ~sgn c\"u/ —--f y 2|g| 

is normalizable. We conclude that 



(2) / \ 
sen c V / 



I 

'vSgn c. 



(6.16b) 



(6.17a) 



(6.17b) 



is a Majorana state with the eig envalue sgnc under the 
particle-hole transformation (5.4). The uniqueness of this 
Majorana state, up to a normalization factor, can be 
proved along the same lines as is done in Appendix ^ □ 



B. The Fermi limit 

The Fermi limit is defined by the condition 

%,D - 0. (6.18) 

In this limit and when the triplet-pair potential van- 
ishes, isolated vortices support finite-energy Caroli-dc- 
Gennes-Matricon bound states in the weak-coupling limit 
Ag//i ^ 1. No Majorana fermions are to be found tied 
to the core of an isolated vortex. 



In the limit ( 6.18| ) with a vanishing singlet-pair poten- 
tial, 



fe{k,k) 

e{k,k) 2^{^'^} 

H.c. -£(fc, k) 

\ H.c. 









-e{k,~k) ) 

(6.19) 

decomposes into the direct sum of the 2x2 Hermitian 
operators 



and 



7^(2) 



e(fc, fc) 



A 



2L 



At 



{0,/c}t -e(fc,/e) 



e{k, k) 
^{z,fc}t -e{k~k) 



.20a) 



(6.20b) 
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We claim that Hamiltonian H^^^ supports two normal- 
ized zero modes if and only if (iff) the chemical potential 
fi>0. 



Proof. Define the short-hand notation 

A, 



(6.21) 



We need 



and 



{z, k} = zk + kz 

= 2kz + [z, k] 
= 2kz + 2i 



{z, ky ~kz + zk 

= 2kz + [z, k] 
= 2kz + 2i. 



(6.22) 



(6.23) 



Equation ( |6.20 ) becomes 
H 

and 

H 



(1) _ f e{k, k) c (kz + i) 
c(fcz-l-i) ~e{k,k) 



(2) ^ ( fc) (^kz _ 
ckz —e(k,k) 



(6.24a) 



(6.24b) 



We are going to show that operator ( |6.24a ) has one and 
only one zero mode iff fi > 0. We will then show that the 
same is true for operator ( 6.24b ). 
We seek a solution to 



U - U^oi- I (1) 



We must solve 



= e(fc,fc)u(i' + c (fc2i9;^ + i) w 



^ c{k2idk+i)u''^'^ ^ e{k,k)v 



(1) 
,(1) 



(6.25) 



(6.26a) 
(6.26b) 



If we take the complex conjugate of Eq. ( |6.26bD , we get 



= e{k,k)u^^> + c{k2id-^ + i} 



,(1) 



(6.27a) 



c(fc2ia, 



7(1) 



e{k,k)v'-^l (6.27b) 



We thus infer that a solution to Eq. ( |6.25| ), if it exists, is 
given by 



,(1) 



±u 



(1) 



where w^^^ is the solution to 



{k, k)u'-l'^ + (±)c {k2id-^ +i)u 



(1) 
± • 



(6.28a) 



(6.28b) 



Zero modes, if they exist, can be labeled by their an- 
gular momentum. We seek a zero mode with vanishing 
angular momentum, i.e., independent of (p. We must 
then solve 



where Mj!^ is the solution to 



< K < oo. 







r(^)u£> + (±)(ic) iKd, + l)u'-i\ 



.29a) 



(6.29b) 



With the help of a global gauge transformation, we can 
always choose c so that 

c := ic (6.30) 

is real valued. Hence, 

e(K)u£^ (k) + (±)5 [kO^ + 1) (k) (6.31) 

with < K < oo admits a purely real or a purely imag- 
inary solution since all coefficients of this first-order dif- 
ferential equation are real valued. We choose the real- 



valued solution. We divide Eq. (6.31) by (±)ck to obtain 
the condition 



0= 9.- + 



l + (±) 



e(K) 



(6.32a) 



whose formal solution is given by 



(1)/ ^ (1)/ \ 

U j_ (K) ~ U j_ [Kq) X CXp 



dK' 



l + (±)' 



(6.32b) 

The formal solution ( 6.32b ) is admissible iff it is nor- 



malizable, i.e., if 



oo 

/ 



< oo. 



(6.33) 



Define 



dn' \^ ^ (±)£i^ 



3.34) 



1_(±)^1 ln^ + (±)^!L^. 
c J Kg 4mc 



For large k. 



so that normalizability imposes the choice 
± = sgn c 



(6.35) 



3.36) 
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and the formal solution (6.32b) becomes 

II. , 



sgn c V / sgn c 

For small k 



<:^(-) = <L(-o)x (f 1" xe-^. (6.37) 



(6.38) 

so that normalizability demands convergence of 





1.6. 



(6.40) 



We conclude that for any choice of gauge such that c = ic 
is real valued 



fsgn 5) ui]li:{K), 



where 



,(1) 



sgn 



< K < OO, 

(6.41a) 



^2 

X e"^^^ (6.41b) 



is a normalizable Majorana mode iff ^ > 0. Observe that 
$gg]^g(fe) is an eigenstate with the eigenvalue — sgnc of 



the particle-hole transformation defined in Eq. (5.4), i.e., 



3.42) 



To show that solution ( 6.41 1 ) is unique, up to a nor- 
malization, one expands Eq. (|6.25 ) in polar harmonics 
labeled by the angular quantum number n G Z (see Ap- 
pendix H). Modes with angular quantum number ±n 
are pairw ise coupled. A formal zero mode of the form 
Eq. (6.41) whereby the function u is substituted by a 
doublet, i.e.. 



U«„(«)=exp -W'^'Jn) V'^^M 



(1) 



t(i) 



with 



dK' 



-(1) 



(6.43a) 



(6.43b) 



and g2^''„(k') both 2x2 matrices, follows. However, it is 
not normalizable. 

It is time to seek a solution to 



= ^(2) 



(2) 



(6.44) 



We must solve 

= e(fc, k)u^^'> + ck2idj^v^'^^ , (6.45a) 
= ck2idkU^^^ ~e{k, fc)^^^) . (6.45b) 

If we take the complex conjugate of Eq. ( |6.45b| ), we get 

0==£(fc,A:)u(2)+cA:2i9^w(2), (6.46a) 
= -cfc2ia;^M(2) - e{k, k)v^'^^ . (6.46b) 



We thus infer that a solution to Eq. ( |6.44 ), if it exists, is 
given by 



(2) 



,(2) 



±U 



(2) 



(6.47a) 



(2) ■ 

where u\_ is the solution to 



.47b) 



When expanding the solution in angular momentum 
modes exp(in(y9), n € Z, the mode n = 1 turns out to 
be the only mode that does not couple to other modes 
via Eq. (|6.47bD. The ansatz 



(6.48a) 



casts Eq. ( 6.47bD in the same form as Eq. (6.31) so that 



(2) 

V 

sgn c 



(6.48b) 



(2) (2) 

Observe that ^loici^) with (fc) given in Eq. 



IS 



an eigenstate with the eigenvalue — sgn c of the particle- 
hole transformation defined in Eq. (^.4[), i.e.. 



(6.49) 



It remains to verify that the spinor ( |6.48 ) is single val- 
ued in re al spac e, i.e., that the Fourier transform of the 
function (B.48a) vanishes at the origin of the complex- 
z = rexp(i0) plane. Hence, we need the small r expan- 
sion of the Fourier transform 

OO 27T 





OO 27T 

= e'' J dnn J dcbe'"^^'""^ x e'"^ wif^ . 



(6.50) 

The K integral is well behaved for large k because of the 
Gaussian factor. Moreover, an upper cutoff to this in- 
tegral can be used up to Gaussian accuracy. If so, a 
Taylor expansion of exp(irK cos (/)) in the integrand can 
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be performed to capture the leading dependence on r. 
The integral over (j) eliminates the term independent of r 
so that 



sgn c 



{r,e) - re'" + 0(r2) 



(6.51) 



is single valued at the origin r ~ Q and thus an admissible 
solution. 

We conclude that there are two Majorana states 






sgnce" 




"sgnc(K) 



(6.52) 

{A and B are real numbers) bound to the core of an 
isolated vortex satisfying the linear profile (^|^) in the 
triplet-pair potential. The Majorana state weighted 
by the coefficient A is related to the Majorana state 
weighted by the coeffi cien t B through the helical sym- 
metry defined in Eq. (5.5) a nd not by the operation of 
time reversal defined in Eq. (5^). This is expected since 
TRS is broken by the vortex. □ 

Lu and Yip in Ref. ^ (see also Sato and Fujimoto in 
Ref. ^5|) also found two Majorana fermions bound to the 
core of a vortex with unit vorticity in a weakly coupled 
(i.e., a large chemical potential compared to the pairing 
potentials) 2D TRS noncentrosymmetric superconduc- 
tor with dominant triplet-pair potential. Their firs t zer o 
mode is the real-space counterpart to the mode (3.41). 
Their second zero mode carries angular moment um n ~ 1 
and is the real-space counterpart to the mode ( ^.48 ). 

In Ref. Qi et al. have studied zero modes bound 
to the core of vortices in TRS ± vpy superconduc- 
tors as well. Viewing the system as a combination of 
a.p^ + ipy superfluid (which corresponds to T-Lmv) and its 
time-reversed partner, a — ipy superfluid (which corre- 

sponds to TiXm), they simultaneously introduced a vortex 
in the former and an antivortex in the latter. In contrast 
to our study of a TRS-breaking vortex, this configura- 
tion of a pair of vortex and antivortex is TRS and the 
two Majorana modes obtained by Qi et al. are connected 
by the operation of time reversal. Whereas the Majorana 
fermions (6.52) arc not robust to a generic perturbation 
that breaks translation invariance, the TRS-protected 
pair of Majorana fermions obtained by Qi et al. is ro- 
bust to any weak perturbation that preserves TRS. 



C. Away from the Rashba-Dirac and Fermi limits 

Majorana fermions tied to vortices are robust to con- 
tinuous changes in the BdG Hamiltonian as long as the 
spectral gap does not close, for all nonvanishing energy 
eigenvalues occur pairwise with the energies ±i? owing 
to the particle-hole symmetry (5.4). There will be one 



(Fig. H). By the same reasoning, region III of Fig. ^ 
does not adujui Majorana fermions bound to the core of 
unit vortices£3 Regions II and III in Fig. ^ differ by the 
even number of Majorana fermions that TRS-breaking 
vortices can accommodate. This distinction is not ro- 
bust to any generic perturbation that breaks translation 
symmetry, for it is not protected by TRS. 



VII. 



DISCUSSION 



(two) Majorana fermion(s) tied to the core of an isolated 
vortex carrying vorticity one in regions I (II) of Fig. ^ 



In this paper, we studied the possible superconducting 
phases of the surface states of 3D TRS topological in- 
sulators and 2D TRS noncentrosymmetric metals. Both 
systems share remarkable magneto-electric effects, how- 
ever their bulk superconducting phases differ in impor- 
tant ways. The difference stems from the topology of 
the bands. Surface states of 3D TRS topological insu- 
lators are topologically equivalent to a single species of 
Rashba-Dirac fermions while noncentrosymmetric metals 
are Fermi like with two Fermi surfaces for large chemical 
potentials. As a result, we find that there is a unique 
superconducting phase in the case of the Rashba-Dirac 
limit while there are two phases in the Fermi limit. 

We studied the phase diagram as a function of the 
strengths of the mean-field pair potentials Ag (singlet) 
and A^. (triplet), as well as ^ (chemical potential). In the 
Rashba-Dirac limit, a single Majorana fermion is bound 
to the core of an isolated and TRS-breaking vortex with 
unit winding number in the superconducting order pa- 
rameter everywhere in the phase diagram in the Aj/Ag-/i 
plane, with the exceptions where the gap closes. The 
gap-closing lines do not separate distinct phases in the 
Rashba-Dirac limit, because one can always connect two 
sides of a gap-closing line by, instead of crossing the line 
directly, going through the point at infinity (Ag = 0) 
without closing the gap. Evidently, gap closing is a nec- 
essary but not sufficient condition to have two distinct 
phases. 

In the Fermi limit, we find instead that there are 
two superconducting phases. These phases correspond 
to singlet or triplet dominated physics. In the singlet- 
dominated phase, we find that an isolated TRS-breaking 
vortex with unit winding number (a full vortex) does not 
bind Majorana fermions. In the triplet-dominated phase, 
we find a pair of Majorana fermions bound to an isolated 
full vortex. Hence, these Majorana states have a distinct 
origin from those found for half vortices in the ± ipy 
superconductors. The physical reason for this difference 
is that TRS forces the spin-resolved pairing amplitudes 
A^^ and A^j^ to be related, and thus one cannot intro- 
duce vorticity in one but not the other, as can be done 
with half vortices in the p^ ± ipy superconductors. 
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Appendix A: Proof of Eqs. (4.4), (4.9), and (4.13) 



All reduced interaction Hamiltonians (4.3)j_(4.8), and (4.12) have summands which can be represented in terms of 
the helicity basis using the transformation ( 2.4d ). To do this explicitly, let ^ and ly run from to 3 with ctq the 2x2 
unit matrix and write 



c c 



^^^^^ S i2^^S fciSj^ I 



(Al) 



The • • • stands for two-fermion contributions that result from anticommuting the operators. 
We first evaluate the matrix products 



cxp 1 



= ± exp 1— ^- 



T cxp ^1 
cxp i 











2 




<^p 








2 




<^p 




^k 




2 




^p 




fk 


2 



COS 







Vfc 




2 




vp 


+ 


fk 




2 




vp 


+ 


Vk 




2 




^p 







ism 



.(1) 



^p 




Vfc 




2 




vp 


+ 


Vfc 




2 




vp 


+ 


"Pk 




2 




^p 




Vk 


2 



(A2) 



Here, we observe that any of the right-hand sides in Eq. (A2) involves one diagonal and one off-diagonal Pauli matrix. 
At this point we can partly settle our constraint to have only Cooper pairs made of electrons of th e same helicity. For 
this case, only products of two diagonal or two off-diagonal terms contribute in the product ( [Al[ ). 
Second, we introduce the notation 



VkxlpX' 



fcA"-fcA'^-pA'"pA'j 



A,A = ±, 



(A3) 



in terms of which we find 



i^i<^^"''4'pi^lk<^^"^- 


-p ^ 




-p 


= Vki 


-|p+ " 


^ Vk- 


-lp- ^ ' 


(A4a) 


i^U^'^pi^lk^^'H^- 


-p ^ 


-i'W^^pi'lk'^^'^- 


-p 


= Vk^ 


-|p- " 


^Vk- 


-lp+ + ' ■ ■ ' 


(A4b) 


4'i'j^°^p4'lk'^^'^- 


-p ^ 


+^y^pi'lk<^^''^^ 


p 


= Vki 


-!p+ " 


-Vk- 


-lp- + ■■■ ' 


(A4c) 




-p ^ 


+i'W^^p^ik^^'^- 


-p 


= — i 




p- 


- Vk-\p+^) ^ • 


(A4d) 



Here, • • • stands for contributions which would lead to Cooper pairs made up of two electrons of different helicity. We 
are left with the task of collecting the phase and trigonometric multiplicative factors from Eq. ( A2 ) . 

For the density-density interaction (O), we have to compute Eq. (Al) with ii = v = 0. According to Eq. (A2) this 
involves collecting the phase and trigonometric multiplicative factors for Eqs. (A4a) and (A4b). We find 



i'^-k^U^-p^-p 



kp 



2 "Pp - i 
COS (%+|p4 



'7fc-|p-J 



Pp-Pk 



,2 ^P 



Vk- 



(A5) 



from which Eq. (4.4) follows. 
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For the Heisenberg interaction (4.8), we have to compute Eq. (Al) with ^ = = 1 , 2, 3. Acco rding to Eq. ( A2) this 
involves coUecting the phase and trigonometric muhiphcative factors for Eqs. ( A4a ) and (A4b). We find 



2 '^P ~ '^fc 



%+|p- + '?fc-|p+ %+|p+ Vk-\p- 



COS 



2 



Vk+lp- + '?fc-|p+ 



(A6) 



sm 



Vk+\p+ + ^fc-|i 



from which Eq. (O) fohows. 

Finally, the Dzyaloshinskii-Moriya intera ction ( 4.12| ) invo lves terms of the type (Al) with (/i, i^) = (2,3), (3,2), 
(1,3), and (3, 1) which in turn lead to Eqs. ( A4c ) and ( A4d ). The calculations yields 



m=l,2 



h,p j.,l—1...3 



= iE^'^"^ {Vk+lP--Vk-\p+)[Di\co, 

k,p 



2 co'^ 2 2 ^"^^ 2 



\Jlk+\p+ - Vk-\p-) 



D'k-p sm 5 sm D^j^^ cos sm 



(A7) 



from which Eq. (4.13) follows. 



Appendix B: Unnormalizability of higher angular momentum zero modes 



W e are going to show that the solutions to Eq. (6.28b) with nonzero angular momentum and the solutions to 
Eq. ( |6.47b| ) with angular momentum other than +1 are not normalizable. We expand 



for Eqs. (6.28b) and ( |6.47b ), respectively. The differential equations mutually couples two and only two angular 
momentum modes. As all coefficients of the differential equation are purely real, the expansion parameters /j."'^„(k) 
and '^^^ be chosen to be either purely real or purely imaginary numbers. Without loss of generality, we make 



(1) 

the former choice. In terms of the doublet Uji „ = 
we find 



f^ni9±.n that represents the two coupled modes labeled by n. 



1 



0) 



(B2a) 



The matrices are given by 



2-{j + n) ±'M 



(B2b) 



Hence, the doublet solution can be written as 



u2i(«)=cxp[-F2i(«)lu2i('^o) 



(B3a) 



with 



(B3b) 
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The demand of normalizability reads 



V%{n)YlJ^£jK) = f dKK [u2)„(^o)]''exp ( - 2¥%{K)Yl]n{K,) < oo. 



(B4) 



In the limit of large k, both matrices g'^^^^k) (j = 1, 2) obey the same behavior. The matrix in the exponent becomes 
for both cases j = 1,2, 



Upon exponentiation it, the condition (B4) then reads 



(B5) 



cosh ( 



(/E(«o)) + (ff2?„(«o)) J T2sinh(^^j/W„(..o)g2;U«o) 
which gives a condition for the initial values 

In the opposite limit of small k, we find for the matrix that has to be exponentiated 



< oo, 



K 



Upon exponentiating it, the condition (B4) reads 



-2 1 



> / d. 



(n + j-l)T2|/2;lM52?«M+0 



n+J-l)^ + - 



(B6) 



(B7) 



(B8) 



(B9) 



All terms that are given explicitly in the curly bracket have to vanish in order to achieve normalizability. This amounts 
to 



(n + j-l)2 + ^-(,i + j-l) 



52?„(«o)- 



(BIO) 



Both condit ions (B7) and ( BIOD are only satisfied simultaneously if n = 1 — j. F or th is mode not to be vanishing, the 
sign in Eq. (B7) must be chosen ± = sgnc, which is only compatible with Eq. (BIO) for ^ > 0. This corresponds to 
the solutions discussed in Sec. ^ and we conclude that these are the only normalizable zero modes for each of the 
blocks 7 = 1,2. 
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